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Abstract. It was known that by a duality transformation, interacting bosons at 
filling factor / = p/q hopping on a lattice can be mapped to interacting vortices 
hopping on the dual lattice subject to a fluctuating dual " magnetic field" whose 
average strength through a dual plaquette is equal to the boson density f — p/q. So 
the kinetic term of the vortices is the same as the Hofstadter problem of electrons 
moving in a lattice in the presence of f — p/q flux per plaquette. Motivated by this 
mapping, we study the Hofstadter bands of vortices hopping in the presence of magnetic 
flux f — p/q per plaquette on 5 most common bipartite and frustrated lattices namely 
square, honeycomb, triangular, dice and Kagome lattices. We count the total number 
of bands, determine the number of minima and their locations in the lowest band. 
We also numerically calculate the bandwidths of the lowest Hofstadter bands in these 
lattices that directly measure the mobility of the dual vortices. The less mobil the dual 
vortices are, the more likely in a superfluid state the bosons are. We find that except the 
Kagome lattice at odd q, they all satisfy the exponential decay law W — Ae~ cq even at 
the smallest q. At given q, the bandwidth W decreases in the order of Triangle, Square 
and Honeycomb lattice. This indicates that the domain of the superfluid state of the 
original bosons increases in the order of the corresponding direct lattices: Honeycome, 
Square and Triangular. When q = 2, we find that the the lowest Hofstadter band is 
completely flat for both Kagome and dice lattices. There is a gap on Kagome lattice, 
but no gap on dice lattice. This indicates that the boson ground state at half filling 
with nearest neighbor hopping on Kagome lattice is always a superfluid state. The 
superfluid state remains stable slightly away from the half filling. Our results show 
that the behaviours of bosons at or near half filling on Kagome lattice are quite distinct 
from those in square, honeycomb and triangular lattices studied previously. 
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1. Introduction 

The Extended Boson Hubbard model with various kinds of interactions, at various kinds 
lattices ( bipartite or frustrated ) at various kinds of filling factors ( commensurate 
/ = p/q or in-commensurate ) is described by the following Hamiltonian [EllfIl[I|: 

H = - t ( b l b i + h - c -) -/^E"* + ^Yl n i( n i ~ !) 

<ij> i i 

+ Vl Y n i n j + V 2 Y U i n k H (!) 

<ij> «ik» 

where n, = b\bi is the boson density and U,V\,V2 are onsite, nearest neighbor (nn) 
and next nearest neighbor (nnn) interactions between the bosons. The • ■ ■ may include 
further neighbor interactions and possible ring-exchange interactions. For a bipartite 
lattice, the sign of t can be changed by changing the sign of hi in one of the two 
sublattices. But in a frustrated lattice, the sign of t makes a difference. 

It is very important to extend Boson Hubbard model in bipartite lattices to 
frustrated lattices such as triangular, dice and Kagome lattices, because of the following 
motivations: 

(1) For atatoms adsorptions on bare graphite, the preferred adsorption sites form a 
triangular lattice. The phase diagrams of coverage ( the filling factor ) verse temperature 
resulting from the competitions of these energy scales are very diverse and rich PQI21- 
It was believed that Eqn^ may capture the main physics of the phenomena. 

(2) Atomic physicists are trying to construct an effective two dimensional frustrated 
optical lattices using laser beams and then load either ultra-cold fermion or boson atoms 
at different filling factors on the lattice. They may tune the parameters to realize 
different phases by going through quantum phase transitions [HUE]. 

(3) In the hard-core limit U — > oo, due to the exact mapping between the boson 
operator and the spin s = 1/2 operator: b\ = Sf ,bi = S~,rii = Sf + 1/2, the boson 
model Eqn^ can be mapped to an anisotropic S — 1/2 quantum spin model in an 
external magnetic field 

H — - 2t Y ( S ? S j + S i S D + Vi E S i S j 

<ij> <ij> 

+ V 2 £ s *s z k -hJ2s? + --- (2) 

«ik» i 

where h — fi — 2V\ — 2V 2 for a square lattice. Note that in this Hamiltonian, 
there is a ferromagnetic coupling in the XY spin components and anti-ferromagnetic 
coupling in the Z spin component. Again, in a bipartite lattice, the sign of t can be 
changed by changing the sign of S*,Sf in one of the two sublattices, but keeping Sf 
untouched, so Eqn. EJ is the same as Quantum Heisenberg Antiferromagnet (QHA). 
However, in a frustrated lattice, the sign of t makes a difference, so Eqn. El is quite 
different from the QHA. The one to one correspondence between physical quantities 
in boson model and those in spin model are the boson density corresponds to the 
magnetization n <-> M, the chemical potential corresponds to the magnetic field \i <r-* h, 
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the compressibility corresponds to the susceptibility k = |^ «-> x = The boson 
number conservation corresponds to the U(l) rotation around z axis, the superfluid 
state < bi >t^ corresponds to the XY ordered state < >^ 0, the charge ordered 
state corresponds to the modulation of < S z >. The supersolid corresponds to the 
simultaneous < >^ and the modulation of < Sf > [§J Ej. In the hard-core 
limit, the EqnQ at half filling ( q — 2 ) has the Particle- Hole ( P-H ) symmetry 
bi <->• b\,Tii — > 1 — rii, it can be mapped to EqnJ2] in zero magnetic field h = with 
the Time-reversal symmetry — > —S^,S- — > —Sf. EqnJT] on triangular lattice at 
q = 2 is the prototype model to study supersolid state with P-H symmetry jjj. 

The model EqnHwith only the onsite interaction on square lattice was first studied 
in Ref. 0. The effects of long range Coulomb interactions on the transition was studied 
in jHj. Very recently, the most general cases in square lattice at generic commensurate 
filling factors / = p/q ( p,q are relative prime numbers ) were systematically studied in 
j^j. After performing the charge- vortex duality transformation the authors in jU] 
obtained a dual theory of Eqn^ in term of the interacting vortices tf) a hopping on the 
dual lattice subject to a fluctuating dual " magnetic field". The average strength of the 
dual " magnetic field " through a dual plaquette is equal to the boson density / = p/q. 
This is similar to the Hofstadter problem of electrons moving in a crystal lattice in the 
presence of a magnetic field [10J. The magnetic space group (MSG) in the presence 
of this dual magnetic field dictates that there are at least g-fold degenerate minima in 
the mean field energy spectrum. The q minima can be labeled as ipi, I — 0, 1, • • • , q — 1 
which forms a q dimensional representation of the MSG. In the continuum limit, the 
final effective theory describing the superconductor to the insulator transition in terms 
of these q order parameters should be invariant under this MSG. If < ipi >= for every 
I = ±, the system is in the superfluid state. If < ipi >^ for at least one I, the system 
is in the insulating state. In the supersolid state one condenses a vortex- 

antivortex pair, but still keeps < ipi >= for every I. In the insulating or supersolid 
state, there must exist some kinds of charge density wave (CDW) ( we assume that every 
boson carries one internal charge ) or valence bond solid ( VBS) states which may be 
stabilized by longer range interactions or possible ring exchange interactions included in 
EqnlTJ Very recently, the dual method was used to study the Extended Boson Hubbard 
model on a triangular lattice [TT] . 

In a recent paper [T2] , one of the authors applied the dual approach of the extended 
boson Hubbard model Eqn^to study the reentrant "superfluid" in a narrow region of 
coverages in the second layer of A He adsorbed on graphite detected by Crowell and 
Reppy's torsional oscillator experiment in 1993 I16j. He showed that there are two 
consecutive transitions at zero temperature driven by the coverage : a Commensurate- 
Charge Density Wave (CDW) at half filling to a narrow window of supersolid, then to an 
Incommensurate-CDW. In the Ising limit, the supersolid is a CDW supersolid; whereas 
in the easy-plane limit, it is a valence bond supersolid. Both transitions are second order 
transition with exact critical exponents z = 2, v — 1/2, rj = 0. The results concluded 
that 4 He lattice supersolid was already observed in 1993. He also applied the same 
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dual method to study H 2 /Kt /graphite system investigated in the recent experiment 
14 and proposed that a judicious choice of substrate could also lead to an occurrence 
of hydrogen lattice supersolid. Implications to the realizations of a lattice supersolid of 
ultra-cold atoms in optical lattices were also given in Ref. [T2"] . 

Note that in the dual vortex picture, there are always interactions between vortices. 
Because the phase factors from the dual magnetic field only appear in the kinetic term, 
the interactions always commute with any generators in the MSG, so will not change 
the symmetry of the MSG. 

In this paper, we study the Hofstadter bands of vortices hopping in the presence 
of dual magnetic field / = p/q on the 5 most common bipartite and frustrated lattices 
such as square, honeycomb, triangular, dice and Kagome lattices. We especially study 
the bandwidth of the lowest bands. There are at least two motivations to study the 
bandwidth of the lowest bands (1) As pointed out in [Oj, as q becomes too large, the dual 
vortex method suffer the following two drawbacks (a) As dictated by the MSG, there are 
q minima in the BZ, so the distance in momentum space between these minima scales as 
1/q, the continuum theory only works at k <C 1/q, therefore applies only at distance ^> q. 
The validity regime of the dual vortex theory shrinks, (b) When integrating out the 
vortex modes away from the minima, one encounters energy denominators determined by 
this bandwidth, so the dual vortex method may completely break down if the bandwidth 
becomes too small. By a simple argument, they estimated that at large q, the bandwidth 
W of the lowest Hofstadter band scales as W ~ e~ cq with c at the order of 1 18] . So the 
smaller the bandwidth, the smaller the valid regime of the dual vortex approach. (2) In 
the dual vortex picture, there are both a kinetic energy term and interactions between 
vortices. The kinetic term favor the moving of the vortices, while the interactions favor 
the localization of the vortices, the competition of the two energy scales may result all 
kinds of phases such as superfluid, CDW, VBS and even supersolid phase [12] . In this 
paper, we focus on the kinetic term only. Calculating the bandwidth of the kinetic term 
is very important, because the smaller of the bandwidth, the more inert the vortices ( 
the less mobil the vortices are ), therefore the boson superfluid state is more likely to 
occur. 

By choosing suitable gauges and solve corresponding Harper's equations in the 5 
lattices, we count the number of bands, determine the number of minima and their 
locations in the lowest Hofstadter bands. The results are listed in Table 1. We also 
numerically calculate the bandwidths of the lowest bands in these lattices at any q and 
test against the estimate W ~ Ae~ cq . We believe that although the argument in 
seems reasonable, it is far from being convincing. So it is important to test this argument 
by quantitative numerical calculations. We find that except the Kagome lattice at odd 
q, the exponential law is indeed satisfied and determine (A, c) for the 5 different lattices. 
The results are listed in Table 2. We find that at given q, the bandwidth W decreases in 
the order of Triangle, Square and Honeycomb lattice. The corresponding direct lattices 
are honeycomb, square and triangular lattices, so the tendency to form a superfluid 
state increases. As shown in the Table 2, when q = 2, the lowest bands in both Dice 
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and Kagome lattices are fiat. In dice lattice, the gap between the second fiat band and 
the lowest flat band is ^/6. It indicates that for the original boson at half filling with 
nearest neighbor hopping on the Kagome lattice, there could be only superfluid state. 
However, in Kagome lattice, the gap between the second dispersive band and the lowest 
flat band vanishes at k — (0, 0), so the second dispersive band can not be ignored even 
in the lowest energy limit. Due to the gap vanishing on the Kagome lattice, we can not 
say definite things about the ground state in the original boson on a Dice lattice. There 
are some previous results on the energy spectra on square, honeycome and triangular 
lattices (THl M U21 EEH UZj with different focuses. Our results on dual Dice and Kagome 
lattices, especially the discussions on the possible boson ground states on corresponding 
direct lattices are new and most interesting. 

There are two equivalent methods to be used to study the Hofstadter bands. One 
is the Magnetic Brillouin Zone (MBZ) method to be employed in the main text. This 
method is physically more transparent and intuitive. Another is the symmetric method 
used in jH] and to be used in the appendix. This method treat x and y coordinates 
on equal footing, so is more symmetric than the first one. In the main text, we will 
use the first method to derive Harper's equations in the 5 lattices and then solve the 
equations analytically at small q and numerically at large q. In the appendix, we will 
use the second method to repeat the calculations. Although the coefficients of Harper's 
equations in the two schemes are different, as expected, we find that they result in the 
same energy spectra. 

In the following, we will first study two bipartite lattices, namely square and 
honeycomb lattices, then we will investigate 3 frustrated lattices namely triangular, 
Dice and Kagome lattices. In the final section, we summarize our results in Table 1 
and Table 2, we also comment on the results on CDW formations in high temperature 
superconductors claimed in P where q as large as 8, 16, 32 are used. In most of the 
cases, we focus on p = 1 case. 

2. Square lattice 

We are looking at the Hofstadter band of vortices hopping around square lattice in the 
presence of magnetic flux / = p/q per square [§| ( Fig.l). In the MBZ method, one 
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Fig 1: (a) a magnetic unit cell of square lattice, (b) Phase factors on bonds, phase factors are not 
shown. 
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magnetic unit cell is q times larger than the conventional unit cell (Fig. la). For the 
simplest gauge chosen in Fig. 1(b), the hopping Hamiltonian is: 

H = -tJ2[\x + ai >< A + \x + a 2 > e i2nfai < x | + h.c] (3) 

X 

In the following, for simplicity, we set t = 1. The eigenvalue equation Hip(k) = 
E(k)if)(k) leads to the Harper's equation: 

- e-^^ik) - 2cos(2vr// + k y )^(k) - e ik ^ l+1 {k) = E(k)^(k) (4) 

where I — 0, ■ • •, q — 1; ^ < k x < -; — n < k y < n. 

For small values of q, Eqn0] can be solved analytically. For large values of q, we 
solve it numerically. There are always q bands. We focus on the lowest energy band 
and its bandwidth. As shown in j^j, there are q minima at (0, 2irfl), I — 0, • • •, q — 1. 
The spectrum for q = 4 is shown in Fig. 2. In order to see clearly all the 4 MSG 
related minima, only part of the energy band close to the 4 minima in the lowest 
band is included. We also numerically calculated the bandwidth of the lowest band 




Fig 2: The lowest energy band of square lattice at q = 4. 



upto q — 18. We found that it indeed satisfy the exponential law W = Ae cq with 
A = 26.05, c = —1.20. In a semi-log plot, it is a straight line which is shown in Fig.3 

What is surprising is that even for the smallest q = 1 which is the no magnetic field 
case, the exponential law is still satisfied. 

3. Honeycomb lattice 

Honeycomb lattice is not a Bravais lattice, it can be thought as a underlying 
parallelogram Bravais lattice with two primitive vectors a\ = x, a 2 = -^x + ^fy plus 
a two point basis located at x + 5 and x + 25 where 5 = |(ax + a 2 ) ( Fig. 2 ). Its 
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Fig 3: The bandwidth of lowest band in square lattice vs q 




Fig 4: (a) A magnetic unit cell of honeycomb lattice (b) Phase factors on bonds, phase factors 
are not shown 

In the MBZ method, one magnetic unit cell is q times larger than the conventional 
unit cell (Fig.4a). In one conventional unit cell, there are also two atoms which are 
labeled by two color indices a and h ( Fig.4). We are looking at the Hofstadter 
band of vortices hopping around a honeycomb lattice in the presence of magnetic flux 
/ = p/q per hexagon. For the simplest gauge chosen in Fig. 4(b), the vortex hopping 
Hamiltonian is: 

H — — t ^2[\x + 5 >< x\ + \x + 5 >< x + d[\ 

+ \x + 5 > e~ i2Tfai < x + a 2 \ + h.c] (5) 
The Harper's equation is: 

- (1 + e i{k * +2nfl) )^{k) - e ik »ip? +1 (k) = E{k)^{k) 

- (1 + e-^ +2 ^)^(k) - e-^Vtx® = E (kmk) (6) 
where I — 0, ■ ■ -, q — 1 is the flavor indices and a, b is the color indices, — - < k x < - . 

For small values of q, Eqn|Hl can be solved analytically. When q = 1, there 
is actually no magnetic field, it is just ordinary tight-binding model. There are 
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two bands: ±^3 + 2(cos k x + cos k y + cos(k x + k y )). The lowest energy band is 

— y/3 + 2 (cos k x + cos k y + cos(k x + k y )). There is only one minimum at (0,0). The 
bandwidth is 3. The q = 2 case is especially interesting, because the original boson 
model can be mapped to a quantum s — 1/2 spin model Eqn|2]in a triangular lattice 

at zero field. For q = 2, there are 4 bands E(k) = ±t\J 3 ± \j2A{k) where A(k) = 

3 + (cos 2ki + cos 2k 2 — cos(2fci — 2k 2 )). The lowest subband is E(k) = —tJ 3 + \j2A(k). 
There are 4 minima at ±(7r/6, — tt/6) and ±(7r/6, 5tt/6). The 4 minima transforms to 
each other under the MSG. 

For general q, there are always 2q bands. As shown in , there are two cases 

(1) q is odd, there are q minima at (0, 2%fl), I — 0, • • -, q — 1. (2) q is even, there are 2q 
minima at (^, — ^ + 27r/Z) where a = ±, / = 0, ■ ■ -, g — 1. 

For large values of q, we solve EqnlH] numerically. Just like in square lat- 
tice, we focus on the energy band near the minima in the lowest energy band. 
The q = 3 and q = 4 spectrum are shown in Fig. 5(a), 5(b) respectively. 




Fig 5: The lowest energy bands of honeycomb lattice at (a) q = 3, (b)q = 4 

We also numerically calculated the bandwidths of the lowest band upto q = 18. We 
found that they satisfy the exponential law W = Ae~ cq with A = 11.82, c = —1.66 for 
both q even and odd. In a semi-log plot, it is a straight line which is shown in Fig. 6. 




Fig 6: The bandwidth of honeycomb lattice vs q 
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What is surprising is that even for the smallest q = 1 which is the no magnetic field 
case, the exponential law is still satisfied. 

4. Triangular lattice 



In the previous two sections, we studied two bipartite lattices. In this section, we study 
the simplest frustrated lattice which is the triangular lattice. As said in the introduction, 
the physics in frustrated lattices could be very different from that in bipartite lattices. 




(a) (b) 

Fig 7:Triangular Lattice (a) magnetic unit cell of Triangular lattice, (b) Phase factors on bonds, phase 
factors are not shown. 

We are looking at the Hofstadter band of vortices hopping around a triangular 
lattice in the presence of magnetic flux / = p/q per triangle. For the simplest gauge 
chosen in Fig. 7b, the Hamiltonian is: 

H = - t Y^[\x + al >< x\ + \x + a 2 > e t27T2fai < x\ 

X 

+ \x + at + a 2 > e^ 2 /( ai +5) < £ | + h.c] (7) 

The corresponding Harper's equation is: 

- 2 cos(fcj, + 2tx fl)^(k) - (e~ ik * + e~ i( ^ +fc! ' +2 " /(2 ^ 1)) )^_ 1 (A : ) 

_ ( e <fc- + e *(*«+*»+2*/(2m)))^ I+1 (jfe) = E(k)ipi(k) (8) 

where I — 0, • • •, q — 1. 

From the phase factors on the bond, it is easy to see that when q is even, there are 
only | unit cells in one magnetic unit cell. Since the magnetic unit cell shrink to |, the 
range of k x in the momentum space double its range accordingly. Therefore in Eqn|Hl 
for q is odd, — - < k T < -, while for q even, — — < k T < — . 

In fact, as shown in [0], there are three cases in a triangular lattice: (1) When q is 
odd, there are q bands. There are q minima at (0, 47r//), I — 0, • • •, q — 1 in the lowest 
band. (2) When q is even, there are q/2 bands. There are still two subcases: (2a). 
q = 2n with n odd, there are q minima at (^, + 47r//), a = ± and I = 0, • • -| — 1. 
(2b). q = 2n with n even, there are | minima at (0, 47r//), I = 0, • • •, | — 1. 

For q = 1 which is the no magnetic case, the energy spectrum is E(k) = 
— 2(cosfca; + cos k y + cos(k x + k y )), there is only one minima at (0,0). For q = 2, as 
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shown in |T2j, the spectrum is E(k) = —2(cos k x + cos k y — cos(k x + k y )) . There are 
two minima located at (±| , 

For large q, we solve EqnJH] numerically. The results for q = 3 ( odd case ), q = 6 
( 2n with n odd case ) and q = 8 ( 2n with n even case ) are shown in Fig. 8a,b,c 
respectively. 




-3.39 



-3.4 
-3.41 



'"TOT" 



(a) 



(b) 




Fig 8: The lowest energy bands of triangular lattice at (a) q — 3, (b)q = 6, (c)q = 8 



We also numerically calculated the bandwidth of the lowest band upto q = 25 for 
q is odd and upto q = 30 for q is even. For q odd, we find A = 9.21, c = 0.82 ( Fig. 9a). 
For q = 2n, both n is odd and even, the bandwidth satisfy the same exponential law 
with A = 55.70, c = 0.83 ( Fig.9b ). 




Fig 9: The bandwidth of triangular lattice vs q (a) q is odd,(b) q is even 
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5. Dice lattice 

The dice lattice is the dual lattice of the Kagome lattice. It can be thought of consisting 
of two interpenetrating honeycomb lattice. Obviously, the dice lattice is not a Bravais 
lattice, it can be thought as a underlying parallelogram Bravais lattice with two primitive 
lattice vectors a\ = x, 0,2 = \x + plus a three point basis labeled a, b, c located at 

— » — * — * 

x, x + 5, x + 25 where 5 = |(ai + a 2 ) ( Fig. 10). In contrast to the honeycomb lattice 
shown in Fig.4, the dice lattice is not a bipartite lattice and has a 3-sublattice structure. 




a i 

(a) (b) 
Fig 10:Dice Lattice (a) magnetic unit cell of Dice lattice, (b) Phase factors on bonds, phase factors 

are not shown. 

We are looking at the Hofstadter band of vortices hopping around a dice lattice in 
the presence of magnetic flux / = p/q per parallelogram. For the simplest gauge chosen 
in Fig 10b, the Hamiltonian is: 

H = -tY^[\x + 6 >< x\ + \x + 5 > e ~ i2w3fai < x + a 2 \ + \x + 5 >< x + d[\ 

X 

+ \x + 26 > e -^ 3 /K+D < x + a 2 \ + \x + 2~S >< x + d[\ 

+ \x + 2~5 > e - i2 -3/( a i+f) < £ + a- + (T 2 \ + h.c] (9) 

The corresponding Harper's equation is 

- (1 + e- t(k » +27r3fl) )^(k) - e'^^ik) - e- i{k y +27Tm+ ^^(k) 

- {e~ ikx + e -i(k x +k y +2*3f{i+i))^c_^ = E(k)ijji(k); 

_ (1 + e -^ +2 ^)^(k) - e ik *W +1 (k) = E(k)tf(k); 

_ e i(k y +2n3f(l+±))^a^ _ ^ e i(fe+fey+27r3/(/+|)) + e ifc -)^ a +1 (^) 

= E(kU(k) (10) 

where I — 0, 1, • • • , q — 1 is the flavor indices and a, b, c are the 3 color indices. 

For the simplest gauge shown in Fig. 10b, we need to distinguish two general cases: 
q = 3n and q 7^ 3n. For q 7^ 3n, we find out there are still two subcases: q is even and 
q is odd. 
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When q is small, we can solve the Harper's equation analytically. For q = 1 which 
is the no magnetic field case, there are 3 bands: ±y^3 + 2(cos/c a ; + cosk y + cos(k x — k y )) 

and 0. The lowest band is E(k) = — ^3 + 2(cosA; x + cosk y + cos(/c x . — k y )). The 
minimum is at (0,0). For q = 2, all the 6 bands are completely flat. The energies 
are E = 0, — y/E, each with degeneracy 2. For q = 3, there are also 3 bands: 
±^6 + 2A(k x , k y ) and 0. The lowest band is E(k) = —\JQ + 2A(k x , k y ) where 
A(k x , k y ) = cos k x + cos k y + cos(k x — k y ) +cos(k y + + cos(k y — k x + ^) + cos(k x + ^). 
The two minima are at (0, 0) and (— 




Fig 11: The lowest energy bands of dice lattice at (a) q — 9, (b)q — 6, (c)q — 5, (d)q = 4 




O 10 20 30 



<<=) 

Fig 12: The bandwidths of dice lattice vs q (a) q = 3n,(b) q ^ 3n and odd,(c) q ^ 3n and even 

In general,there are four cases in the dice lattcie: (1) q = 3n,^p < k x < — . There 
are g bands. We also need to distinguish two subcases: (la), n is odd, there are 2n 
minina at (0, ^l) and (— |^ + — Z), i = 0, • • -, n — 1. q = 9 case is shown in Fig. 11(a). 
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(lb), n is even, there are n minima at (— |^ + ^l), I = 0, n— 1. g = 6 case is shown 
in Fig. 11(b). For both cases, the bandwidth falls as 14.73e~°' 55q as shown in Fig. 12 (a). 
(2) q ^ ■^> n i~f' < k x < ^. There 3q bands. We also need to distinguish two subcases: 
(2a). q is odd, there are 2q minima at (^p-, — tt + yL + ^i/) a = ±, / = 0, • • •, g — 1. g = 5 
case is shown in Fig. 11(c). The bandwidth falls as 0.54e~ a54lJ as shown in Fig. 12(b). 
(2b). q is even, there are q minima at (|, | + y/), / = 0, • • •, q — 1. q = 4 case is shown 
in Fig. 11(d), the bandwidth falls as 1.0 le -0 549 as shown in Fig. 12c. 

It seems to us that all the four cases have the same c within numerical errors, but 
with different magnitudes A. 

6. Kagome lattice 

Kagome lattice is not a Bravais lattice either, it can be thought as a underlying 
parallelogram Bravais lattice with two primitive lattice vectors c?i = x, a 2 = \x + 
plus a three point basis labeled a, b, c located at x, x + di/2, x + a 2 /2 as shown in Fig.13. 
Note that the Kagome lattice contains both triangles and hexagons. 




Fig 13: Kagome lattice (a) magnetic unit cell of Kagome lattice, (b) Phase factors on bonds, phase 
factors are not shown. 

We are looking at the Hofstadter band of vortices hopping around a Kagome lattice 
in the presence of magnetic flux / = p/q per triangle and 6/ flux quantum per hexagon. 
So overall, there are 8/ flux quanta per parallelogram. For the simplest gauge chosen 
in Fig 13b, the Hamiltonian is: 

H = -t^2[\x+ ^ >< x\ + \x + ^ > e i2w8fai <x\ + \x + d[ ><x + ^-\ 

x 

+ \x + a 2 >< x + — | + \x + — + a 2 > e tZ7Tj < x + on + — | 

+ |f + I > e ^f(a 1+ l) < f + ^| + h.c] (11) 
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The corresponding Harper's equation is: 

- tf(k) - e^^Uik) - (e ik y + e-^ l )^{k) = E(k)^(k); 

- if>f(k) - e- ik *4>? +1 (k) - e- l27r8f(l+ ^^{k) - e-^ 27r/+fe -^Vf+i(^) 
= E(k)^(k); 

- (e~ ik y + e l27r8fl )^(k) - e l2 " 8/( ' + %f (k) 

_ eW+^-ky)^^ = E (k)^(k) (12) 




For the gauge chosen in Fig. 13a, we can solve the spectra at q — 1, 2, 4, 8 exactly, 
because for all the four cases, we only need to solve a 3 by 3 matrix whose 
secular equation is a cubic equation A 3 — (4 + A(k x , k y ))X + 2 cos( ^-)A(k x , k y ) = 
where A(k X) k y ) = 2 + 2(cos/c x + cosk y + cos(k x — k y )). There are 3 bands. For q = 1 
which is the non-magnetic case, the 3 bands are— 1 ± yT + A(k x , k y ), 0. The lowest band 
is E{k) = — 1 — \Jl + A{k x , k y ) whose minimum is at (0, 0). For q = 2, the 3 bands are 
—2, 1 — \Jl + A(k x , k y ), 1 + \Jl + A(k x , k y ). We can see the lowest band is completely 
flat, the second band touches the lowest band at k — (0, 0) where the gap vanishes ! For 
q — 4, the 3 bands are ±^4 + A(k x , k y ), 0. The lowest band is E(k) = —^4 + A(k x , k y ). 
The minimum is at (0,0). For q = 8, we need to solve the cubic equation numerically, 
the minimum of the lowest band is found to be at (0, 0). 

In general, there are 5 cases in Kagome lattice: (1) q — n is odd, — < k x < -, 
there are q minina in the spectrum at (0, —I), I — 0, 1, • • • , q — 1, q — 3 case is shown 
in Fig. 14(a). We find that the bandwidth does not satisfy the exponential law as 
shown in Fig. 15a. (2) q = 2n with n odd, < k x < — , there are | minima at 
(0, — /),/ = 0, 1, • • • , q/2 — 1, q = 6 case is shown in Fig. 14(b). But when q = 2, 
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as shown above, the lowest energy band is completely flat. From Fig. 15b, we can 
clearly see two separate straight lines. We divide the data into separate sets. For 
set 1 in Fig.l5bl, the bandwidth falls as 0.33e -a2 ° 9 . For set 2 in Fig.l5b2, we have the 
bandwidth falls as 0.10e~°' 20lJ .(3) q = An with n odd, < k x < — , there are | minima 
at (0, —I), I = 0, • • -,q/A - 1, q = 12 is shown in Fig.l4(c). The bandwidth falls as 
2AQe~°- 21q as shwon in Fig.l5(c). (4) q = 8n, < k x < I = 0, 1, • • - , q/8 - 1, there 
are also two subcases (4a) when n is even, there are 2n minima at (— f^, — ecu + + —I) 
a — ±. q = 16 is shown in Fig. 14(e). (4b) When n is odd, there are n minima at (0, —I), 
q = 24 is shown in Fig. 14(d). For both cases, the bandwidth falls as 13.87e -0 ' 2113 as shown 
in Fig. 15(d). There are 3n bands in all these cases. 

It seems to us that all the five cases have the same c within numerical errors, but 
with different magnitude A. 
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-12 



-16 



8 16 

q 



40 

q 



(a) 



(b) 



y= -0.20q-1 .12 




(b.1) 




(b.2) (c) (d) 



Fig 15: The bandwidth of the bands in Kagome lattice at (a) q is odd, (b)q — 2n with n odd, 
(b.l)g = 2n case 1, (b.2)q = 2n case 2, (c)q — An with n odd, (d) q = 8n 

7. Summary and conclusions 

In this paper, we have studied the energy spectra of the Hofstadter band of vortices 
hopping on five lattices in the presence of magnetic flux / = p/q per smallest plaquette. 
Our results on dice and Kagome lattices are most new and interesting. The number of 
the energy bands and the number of minima in the lowest band in the five lattice are 
listed in table 1. 
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bipartite 


square 


Honeycomb 




K = v 




2n, n = e 
n, n = o 


KJ l J f ill l. J. o 


2q bands 


frustrated 


triangular 


dice 


Kagome 


N 2n = \n, n = e ^ 

* 1 9n 71 — n 
i —II. it — \j 


N*> = < 


n, n = e 
2n, n — o ' 


jV™ = n = o 


| bands 


3n bands 


Nq n = n = o 


N° = q 


N q^3n = Q , 


Nq U = n = o 


q bands 


N°^ n = 2q 


N sn = | 2n, n = e 
q In, n = o 


above two have 3q bands 


all cases have 3n bands 



Table 1: The number of minimum of the lowest Hofstadter bands in the 5 lattices. N™ means 
q = n. Suffix e and o mean even and odd. We also list the total number of bands just below each cases. 
At q = 2, the lowest band in Kagome and dice lattices is completely flat. 

It was argued in |5j that for large q, the bandwidth of the lowest energy Hofstadter 
band in square lattice scales as W = Ae~ cq with c ~ 1. We believe that although the 
argument seems reasonable, it is far from being convincing. So it is important to test 
this argument by quantitative numerical calculations. We tested the rule by numerically 
calculating the bandwidths of the lowest Hofstadter bands in the 5 lattices. We found 
that this rule is indeed satisfied for all the lattices even for smallest values of q — 1 
except the Kagome lattice for q is odd. For Kogome and dice lattices, the lowest band 
is completely flat at q — 2. The results of (A, c) for the five lattices are listed in the 
table 2: 



Square 


Honeycomb 


Triangle 


Dice 


Kagome 


(26.05, 1.20) 


(11.82,1.66) 


(55.70, 0.83) e 


(14.73, 0.55) 3n 


odd, not apply 






(9.21,0.82) o 


(1.01, 0.54)e^n 


(0.33, 0.20) 2n>1 
(0.10,0.20) 2n>2 ' 








(0.54,0.54) o ^n 


(2.46, 0.21) 4 „ 










(13.87, 0.21) 8n 



Table 2: The bandwidths parameters (^4, c) of the lowest Hofstadter bands in the 5 lattices. Suffix 
e(o) means q is even (odd). At q = 2, the lowest band in Kagome and dice lattices is completely flat. 

From table 2, it is easy to see that for a given lattice, although the prefactor A 
could be different for different cases, the c remains the same within the numerical errors 
for a given lattice. For Kagome lattice when q is odd, the bandwidth does not satisfy 
any exponential decay law. However, in any other cases, they do satisfy the exponential 
laws. The peculiarity of Kagome lattice may be related to the fact that, in contrast to 
all the other 4 lattices, Kagome lattice has both triangles which enclose / flux quanta 
and hexagons which enclose 6/ flux quantum. 

As said in the introduction, the first motivation to study the bandwidth is to look 
at the valid regime of dual vortex approach in the 5 lattices. If q is too large, the 
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bandwidth becomes too small, the dual vortex approach may not be valid anymore. For 
example, on square lattice, when q = 4, W = 0.21 is already very small. This fact puts 
some doubts on the results of CDW formations in high temperature superconductors 
in P where q as large as 8, 16, 32 are used. In fact, large q means very dilute boson 
density in the direct lattice. In this case, the superfluid is probably the only ground 
state anyway except there are very very exotic long range interactions in Eqn^ which 
may stabilize CDW and VBS. Fortunately, the q — 2 ( which is the smallest non-trivial 
case ) in honeycomb lattice was applied by one of the authors to study Helium and 
Hydrogen adsorption problems on various substrates in jTJ]. 

The second motivation is to study the tendency for interacting bosons to form a 
superfluid in the 5 lattices. The bandwidth is proportional to the vortex hopping matrix 
element, so the smallness of bandwidth favors the localization of the vortices, therefore 
enhance the tendency to form a superfluid. At given q, the bandwidth W decreases in 
the order of Triangle, Square and Honeycomb lattice. The corresponding direct lattices 
are honeycomb, square and triangular lattices whose coordination numbers are 3, 4 and 
6. It is known the the higher coordination, the easier for bosons to get the ordered 
superfluid state. The conclusions achieved in dual lattice are indeed consistent with our 
intuition in the direct lattice. 

As shown in the table 1, when q = 2, the lowest bands in both Dice and Kagome 
lattices are flat. In dice lattice, the gap between the second flat band and the lowest 
flat band is \^6. However, in Kagome lattice, the gap between the second dispersive 
band and the lowest flat band vanishes at k = (0, 0), so the second dispersive band can 
not be ignored even in the lowest energy limit. In dice lattice, all the three bands are 
flat, the vortices are completely inert, the interactions certainly favor the localization 
of the vortices. It indicates that for the original boson at half filling ( q — 2 ) with 
nearest neighbor hopping on the Kagome lattice, there could only be a superfluid state. 
Slightly away from half filling, it was known that the superfluid state is stable against 
small number of vacancies or interstitials, we expect the superfluid state remains stable. 
This is in sharp contrast to bosons at half filling hopping on triangular lattice where 
there is a dispersion in the lowest band as shown in this paper. Due to the competition 
between the kinetic energy and the interactions between the vortices, there is a transition 
from a superfluid to a supersolid state as shown in [Zj. The q = 2 case at square |Hj 
and honeycomb jT5j lattices were shown to have CDW or VBS to superfluid transition 
in Ising or easy-plane limit. Slightly away from half filling, in the CDW or VBS side, 
it was shown in J2| that there must be a CDW or VBS supersolid state intervening 
between commensurate CDW or VBS to In-commensurate CDW or VBS in Ising or 
easy-plane limit. Obviously, the behaviors of bosons on Kagome lattice at or near half 
filling ( q = 2 ) are quite distinct from those in square, honey and triangular lattices. 

It is important to stress that the exactly flat bands at q = 2 at Dice lattice are 
completely due to the special lattice structure of Dice lattice which localize the vortices. 
The dual vortex theory immediately leads to the boson superfluid state in the Kagome 
lattice. Of course, the bandwidth goes to zero at large q in any lattices. However, as 
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stated in previous paragraphs, the dual vortex theory is not valid any more at sufficiently 
large q. 

The main body of this paper only discusses the p = 1 case. Taking complex 
conjugate on the Harper's equation Hip(k) = E(k)ip(k) leads to H*ip*(k) = E(k)ip*(k). 
Obviously, H* corresponds to — / which is equivalent to 1 — / = 1 — p/q, so p — 1 has 
the same energy spectra as p = q — 1 for non-interacting vortices. Of course, vortex 
interactions will not have the periodicity / — ► 1 + / anymore, so p = 1 and p = q — 1 in 
Eqn^niay not be equivalent. 

In a future publication, we are trying to understand by analytical methods similar 
to the ones used in [2H] (1) why W = Ae~ cq is satisfied at even smallest value of q 
? Is this a unique feature of any tight binding model ? (2) For different cases on 
triangular and dice lattices and Kagome lattice for q is even listed in Table 2, why c 
is the same within the numerical errors, while A differs ? (3) What is the bandwidth 
rule in Kagome lattice for odd q ? We will also construct MSG's for dice and Kagome 
lattice to understand the energy spectra structure in Table 1. 

J. Ye thanks E. Fradkin for helpful discussions. 

APPENDIX 

In this appendix, we simply list the Harper's equations in the symmetric method first 
used in JH]. They may look different from those corresponding equations got by MBZ 
method used in the main text, but we show that both lead to the same Hofstadter bands 
in the five lattices. This check ensure the correctness of the results in the main text. 

(1) Square lattice 

- e- ife «Vz-i " 2 cos(k x + 2nfl)^ - e^J+i = Eipi 

(A-l) 

where I = 0, • • •, q - 1; =f < k x < J 

(2) Honeycomb lattice 

- (1 + e *fe+ 2 -/0)^ _ e ik y^ +1 = Eift 

- (1 + e ~KkMfi)^b _ e -<*»^,6_ i = E ^a 

(A-2) 

where I = 0, • • •, q - 1; - J < k x < J. 

(3) Triangular lattice 

- 2 cos(k x + 2irfl)ih = Eipi (A-3) 

where I = 0, • • •, q - 1, < k x < s for q odd, I = 0, • • •, q/2 - 1, < k x < 2s for q 
even. 
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(4) Dice lattice 



- (1 + e i(k*+^3fl)^a _ e <fe^a +i = 

_ e i(k x +2n3fl)^a _ ^i(k v -2tt /) + e i(k x +k y +2ir3f(l+i)-4irf) = j^c. 

- (1 + e -i(*«+2*3/I)^6 _ e" ifc ^f_ 1 - e -*(^+2-3/Z)^c 

_ ^-iiky-^f) _|_ e -i(k x +ky+2ir3fl-4wf)j ll i J c ^ _ j^^a (A-4) 



where for g 7^ 3n, Z = 0, • • •, q — 1, < k x < -. For q = 3n, I = 0, ■ • • , q/3 — 

1, < jfe < 3ZL 

(5) Kagome lattice 



- (1 + eiVto+^Mtyf - e - i2 * f 4) C l+1 - e -<(*r/+*y-*«-^8/l)^c = ^6. 

_ (1 + e -i(fc+Jhr8/0^6 _ ^,c +i _ e -<*v^,c = (A-5) 

where for q is odd, Z = 0, • • •, q — 1, — < k x < -, for g = 2n with n odd, / = 
0, ■ ■ -, q/2 - 1, — < k x < —. For An with n odd, I = 0, • • •, q/A - 1, < Jfe* < ^. For 
g = 8n, Z = 0, ■ -\ g/8 - 1, % < ^ < f. 
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